Toward engineered quantum many-body phonon systems 



0. 0. Soykal and Charles Tahan 

Laboratory for Physical Sciences, 8050 Greenmead Dr, College Park, MD 20740 

Arrays of coupled cavity quantum phono dynamical systems in silicon are considered. We study 
physical systems that can exhibit, e.g., Mott insulator states of phonons due to a strong phonon- 
phonon interaction (which is mediated by the impurity-cavity-phonon coupling). Our results indicate 
that quantum many-body phonon systems are achievable both in on-chip nanomechanical systems 
in silicon and DBR phonon cavity heterostructures in silicon-germanium. Experimental procedures 
to detect these states and temperature considerations are given. Cavity-phoniton systems enable 
large phonon-phonon interactions while single phonons offer long wavelengths for forming extended 
quantum states; they may have some advantage in forming truly quantum many-body mechanical 
states as compared to other optomechanical systems. 
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Polaritons in coupled-cavity arrays have received great 
interest for studying strong-correlations and collective 
behavior in light-matter systems pJB]. Simultaneously, 
nanomechanics and optomechanics are driving toward 
the truly quantum regime of mechanical systems [MS], 
where, for example, single photons interact with the low- 
est mechanical mode of a resonator. It is natural to con- 
sider whether these latter systems could exhibit many- 
body quantum interactions in new configurations, allow- 
ing for quantum many-body mechanical systems. To pro- 
vide coherent and strong interaction between mechanical 
modes in a controlled way requires a non-linearity, how- 
ever, analogous to photon blockade [TO] . 

Optomechanical coupling — between photon and me- 
chanical modes — may provide one avenue [8[ |9] to pro- 
duce quantum many-body mechanical systems. At 
present, however, optomechanical coupling must improve 
by a factor of ~ 140 to reach the quantum limit [9j [TT] . 
Also, mechanical resonators, considered as a quantum ob- 
ject, have very short de Broglie wave lengths because of 
their mass, limiting the potential for extended quantum 
states. An alternate system to enable strong coupling 
has been proposed for acoustic phonons p~2j [13] , where a 
cavity phonon hybridizes with a semiconductor two-level 
system (TLS) providing a true analogue to the cavity- 
polariton called a cavity-phoniton, which can easily enter 
the strong coupling regime. In addition, bulk-like single 
phonons in silicon can have long thermal de Broglie wave 
lengths, enabling extended quantum states. 

In this Letter, we introduce two experimentally feasi- 
ble systems where man-made many-body phonon states 
can be realized. We begin by identifying the physical 
parameter regime where many-body Jaynes-Cummings- 
Hubbard Hamiltonians [T4HT6] are realizable, finding 
such phases as, e.g., the Mott insulator states ("Mott 
lobes"). Then, as a starting point for considering real 
experimental setups, we consider a finite array consist- 
ing of only two cavity-TLS sites, calculating the super- 
splitting, the phonon blockade effect, and the response 
to the driving field strength which would be seen in a 




Figure 1. (color online) Schematic of a strain-matched sili- 
con superlattice heterostructure (acoustic DBR with layers of 
Si x Ge(i_ ;r )/Si 2/ Ge(i_ y )) containing donors is shown (left). A 
two-dimensional phononic crystal structure with acceptors 
placed at the cavity sites is also shown (right). 

measurement. We conclude by considering larger system 
sizes, showing that extended arrays behave fundamen- 
tally differently than the small two-site model under the 
same hopping and driving field conditions. 

Schematics of two possible realistic device designs are 
shown in Fig. 1. Our first device proposal involves 
the acoustic phonon cavities constructed from DBR het- 
erostructures via alternating layers of Si^Gei-^ [T7HT9] . 
These structures can be further engineered to possess 
multiple Si cavity regions in a row. In such a setup, 
the overall reflectivity of the layers between any two Si 
cavities simply relates to the phonon inter-cavity hopping 
frequency tij . A suitable donor placed in each of these Si 
cavities can be strongly coupled (a regime where coupling 
frequency is much larger than the donor relaxation and 
cavity loss rates, g ^> T, k) to a specifically chosen single 
cavity-phonon mode u, hence, forming a cavity-phoniton 
[T2] . Our second device design is directly borrowed from 
the concept of nano opto- mechanical phononic crystals. 
An engineered disturbance in periodicity can be used for 
trapping a desired phonon mode in a given region. Place- 
ment of an acceptor impurity into each of these regions 
[T3] will lead to cavity-phonitons with engineered inter- 
cavity tunneling. 



To determine the parameter range for hopping and 
transition frequencies of quantum phase transitions, we 
first consider an equilibrium system where the phoniton 
number density is fixed. This is a good approximation 
where the phoniton lifetime is longer than the thermal- 
ization time. For phoniton arrays consisting of phospho- 
rus donors (or boron acceptors) and phonons in a silicon 
phononic crystal or a DBR array (see Fig. 1), the total 
many-body Hamiltonian is given by the standard Jaynes- 
Cummings-Hubbard (JCH) model [3 EHU H3 EE] , 
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where ai(a\) is the phonon annihilation (creation) oper- 
ator at a given cavity site z, whereas af(a~[) is the exci- 
tation (de-excitation) operator of the donor at that site. 
The inter-cavity phonon tunneling is given by the hop- 
ping frequency for the nearest neighbor cavity sites i 
and j. The regular Jaynes-Cummings Hamiltonian Hjc 
corresponds to the interaction of a single mode of the cav- 
ity phonon with a TLS [22]. The fast oscillating terms 
(i.e. cr^a\) responsible for virtual transitions have been 
dropped via rotating wave approximation. The third 
term in Eq. Q is solely responsible for an effective, non- 
linear on-site phonon repulsion in analogy with photon 
blockade EQI23]. 



The phase transition between a Mott insulator (MI) 
and a superfluid phase (SF) can be determined in the 
grand canonical ensemble where a chemical potential \i 
introduced as 1-L = Hjch — M zCi N{ fixes the number 
density. The operator N = J2i Ni = J2i a l a i + a i a i 
defines the total number of excitations. For simplicity, 
one can assume the random on-site potential with zero 
mean (e.g. fluctuations of the chemical potential), 
vanishes and is assumed to be a uniform short-range 
hopping [3, 20J. The boundary between the MI and the 
SF phases (Mott lobes) is determined by the value of \i 
for which adding or removing a particle does not require 
any energy. Introducing the SF order parameter, ip = 
(a^ via mean- field theory and applying the decoupling 
approximation, i.e. a\dj = (a\)aj + ctl( a j) ~ ( a !)( a j)? 
[T4] we obtain the mean-field Hamiltonian, 
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The correlation number z is the number of nearest neigh- 
bors in a given array geometry. Minimization of the 
ground state energy E of the mean-field Hamiltonian for 
different parameter ranges of /i, cj, and t for phosphorus 
(donor) and boron (acceptor) in silicon yields the Mott 
lobes in Fig. [2^i-b. For the calculation of Mott lobes, in 
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Figure 2. (color online) (a-b) For a many-body phonon-qubit 
system involving P:Si donors (left) and B:Si acceptors (right), 
the SF order parameter, ip, is shown as a function of the 
phonon hopping frequency t and chemical potential \i with 
cavity frequency of uj. MI lobes corresponds to the regions 
of ip = (blue) where number of phonons in each lobe is 
constant ((n) = 0,1,2,...). SF phase corresponds to tp ^ 0. 
(c) Thermal average phonon number per site is shown for 
various temperatures at zero hopping. Plateaus of constant 
(n) corresponds to MI states. 



the case of phosphorus donor impurity, acoustic DBR de- 
sign with correlation number z = 2 (Fig. [I]) is used. The 
donor valley states ls(Ai) and ls(T2) make up the two- 
level system with transition frequency of e = 0.7 THz 
corresponding to a wavelength of roughly A w 12 nm 
[T2] . Due to this small wavelength, DBR heterostruc- 
tures capable of small cavity lengths are the most suit- 
able device structures for maximal coupling. Hence, the 
large array of silicon/DBR heterostructure phonon cav- 
ities can be designed to support a fundamental longi- 
tudinal acoustic (LA) phonon mode in resonance with 
the donor transition (uo = e). In the case of the boron 
acceptor impurity, the transverse acoustic (TA) phonon 
modes of the cavities are reported to yield the maximum 
coupling [13]. TA phonon cavity mode of uj = 14 GHz 
(A = 390 nm) is needed to be in resonance with the spin 
splitting (in the presence of a uniform magnetic field of 
B = 1 T) of the boron valence band acting as a TLS. 
However, at this large wavelength, DBR phonon cavi- 
ties are more difficult to construct due to the critical 
thickness constraint [24], and 2D phononic crystal de- 
signs [6 need to be implemented. For our calculations, 
we used a quality factor of Q = 10 5 currently achievable 
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Table I. Parameters used for a cavity phonon-TLS pair con- 
sisting of a phosphorus (P) donor or a boron (B) acceptor in 
silicon. 



by both designs. The thermal average phonon number 
(n) per site versus \i for various temperatures is shown 
in Fig.^. It is defined by (n) = 1/Z ± ne- En >±/ kBT 

where E Ui ± = (uj - /i)n + (A ± a/A 2 +4^ 2 n)/2 are 
the energy eigenvalues of Hmf with zero hopping and 
Zq = Tr[e~^ MF//kBT ] is the grand canonical partition 
function for the unperturbed (t —> 0) system. The stable 
MI states (compressibility, d(n)/dfi = 0) quickly shrinks 
with increased temperatures. The maximum temper- 
ature allowed to access the first MI state is given as 
T = 0.04-0.06 g/ks in terms of coupling strength. 

The MI and the SF states exhibit different coherence 
characteristics which can be accessed via coherence (cor- 
relation) function measurements [25 , 26 in setups similar 
to modified homodyne/heterodyne or Hanbury-Brown- 
Twiss. Each of these techniques generally require single 
phonon detectors. However, even with single phonon de- 
tectors unavailable, another useful tool, a so called the 
phonon-to-photon translator (PPT), can be deployed to 
coherently convert phonons to photons, therefore allow- 
ing the optical detection techniques to be applied on the 
cavity phonon-TLS system if necessary [T3l [27] . 

The many-phoniton system is driven at each site by a 
phonon field of amplitude Vti and frequency uod- Switch- 
ing to the rotating frame of the driven field yields the 
time-independent Hamiltonian given by 
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where Ae = e — uJd (Acj = uj — UJd) is the detuning be- 
tween the driving field and the TLS (cavity). The driving 
field Hamiltonian is separated as Hg. In the case of dis- 
sipation defined by the cavity loss rate (k) and the qubit 
relaxation rate (T), the master equation for the density 
matrix is given by 
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Figure 3. (color online) (a) Energy schematics of a single, 
two coupled, and an infinite array of coupled cavity phonon- 
donor pairs are shown, (b) For two coupled cavities each 
containing a resonant phosphorus P donor, the transmission 
amplitude versus the detuning between the coherent drive 
(ojd) of strength Q = 2Q C and the cavity field (uj) is shown 
for hopping frequencies of t = 0,t = 0.2g. For a weak 
drive (Q < fi c ) and zero hopping t = 0, system exhibits 
a Lorentzian response (green line) (c) The second-order co- 
herence function 7^ versus the hopping frequency for drive 
strengths of Q — 2Q C and Q — 5Q C , both in resonance with 
the LP branch (ujd — uj — g — t). Donors are detuned by the 
hopping bandwidth A = t and in resonance with the anti- 
symmetric cavity-phonon mode, (d) For an infinite array, 
7^ versus hopping for O < J] c and Q — 2Q C . Donors de- 
tuned by A = 2t. (e) Experimental read-out scheme from 
a single site by a homodyne/heterodyne or modified HBT 
setup. 



where the Lindblad super operator is defined as L[0]p = 
6p6 ] - {O f O, p}/2 PB]. The number of elements of the 
density matrix pij needs to be determined from Eq. ([5| 
are given by (2(A + l)) 2nc , where n c is the number of 
cavities with a single donor /acceptor inside. 

First, we examine the single phoniton system under 
different driving field and hopping conditions. This can 
be done by driving and measuring the heterodyne ampli- 
tude of a single site in the case of zero hopping (t = 0) 
and resonance (s = uj). As seen from the Fig. [3]o (green 
line), for weak driving field strengths smaller than the 



critical value Q < Fl c = (k + T)/4 [29], the system ini- 
tially lies in the linear response regime and it exhibits a 
Lorentzian response to the driving field frequency. The 
critical coherent drive strength is estimated as ~ 42 
MHz and ~ 2 MHz for a phoniton composed of phos- 
phorous donor and boron acceptor, respectively. With 
increasing field strengths, this response breaks down and 
a super-splitting [29] of the phonon field amplitude oc- 
curs (blue line). Intuitively, this behavior can be un- 
derstood as a coupling of the driving field only with 
the antisymmetric dressed state ((|0,e) — |l,g))/\/2) 
and the ground state |0,g), therefore forming a two- 
level system (TLS). TLS treatment will stay valid with 
the driving field strength as long as the non-linearity of 
the Jaynes-Cummings Hamiltonian will only allow single- 
phonon excitations, preventing access to the higher mul- 
tiple excitation manifolds, therefore, causing a phonon- 
blockade. In a single cavity system, the lowest two and 
single excitation energies are given by 62 = 2u — gV2 and 
ei = oJ — g, respectively. This yields to the necessary con- 
dition Qi <C g(2 — \^2) (Qi <C 62 — 2ei) of single excitation 
only subspaces of the system, also known as the "dress- 
ing of the dressed states" [30 , 31 . As the single phoniton 
system still exhibits super-splitting (Q = (tt+r)/2), turn- 
ing on the hopping parameter (t = 0.2g) makes the two 
phoniton states (one phoniton in each cavity) available 
to occupation. This results with a clear shift in eigen 
frequencies and an appearance of a third peak (red line) 
at the heterodyne amplitude spectrum. 

The second order coherence function is defined by 
7< 2 >(0) = (a^a^aa)/(a^a) 2 = (((An) 2 ) - (n)) /(n) 2 + 1, 
where the variance is An = n — (n). MI phase is 
identified by 7 (2) (0) = 1 - l/(n) < 1, whereas SF 
phase by 7^(0) = 1 [32 . For the two coupled phoni- 
ton case, we calculated the second-order coherence func- 
tion 7^) versus the hopping frequency for different field 
strengths. Through out all hopping frequencies, qubits 
were kept detuned from their encapsulating cavity mode 
by A = — e = tto ensure a resonance with the antisym- 
metric mode (lowest) of the overall coupled cavity mode. 
At this detuning choice, the eigen energy difference be- 
tween the ground state (GS) and the lower phoniton (LP) 
branch is given by a simple relation AE = uo — g — t. The 
driving field always kept in resonance with this energy 
difference uod = AE to simulate a TLS system. However, 
for resonant driving purposes, this detuning is not neces- 
sary, as long as one can determine the energy difference 
between GS and LP each time hopping and/or coupling 
parameters are changed. As shown in Fig.[3]3, even in the 
case of strong driving field, ft ^> Q c , the two phoniton 
system exhibits a phonon anti-bunching. 

For large cavity arrays, the mean-field theory and den- 
sity matrix master equation can be applied together for 
weak coherent drive and strong coupling regime [25j [33] . 
Starting from the Hamiltonian in Eq. Q, application of 
the mean field ip = (a) and decoupling approximation 
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in the presence of a coherent driving field. Including 
the cavity loss and qubit relaxation, the master equa- 
tion is same as Eq. ([5| with only driven system Hamil- 
tonian J#s is replaced by the mean-field Hamiltonian 
The SF order parameter ip is evaluated by the 
self-consistency check ip = Tr(pa). For phonitons com- 
posed of P donors, we calculated the second-order coher- 
ence function 7^) versus the hopping frequency for two 
different field strengths, Q <C ^ c and Q = 2Q C = 84 MHz 
in Fig. [3]i. For our particular donor choice, the critical 
drive strength is much smaller than the coupling strength 
^c/g ~ 0.006 due to already small amounts of donor re- 
laxation and cavity loss. For a boron B acceptor, the ra- 
tio is estimated as Q c /g ~ 0.094. The infinite phoniton 
array exhibited a smooth transition from incoherent to 
coherent case, as expected, indicating a phase transition 
from MI to SF state by increasing the hopping frequency. 

We have considered the properties of arrays of quan- 
tum phono-dynamical systems in silicon and show that 
small arrays will demonstrate new behavior and are re- 
alizable and measurable with present techniques. The 
observation of QPTs in large arrays will likely require 
extremely low effective temperatures (at least within 
the approximation considered here [34]), i.e., for P:Si, 
T = 2-3 mK (g = 1 GHz), and for B:Si, T = 40-60 /iK 
(g = 21 MHz). (Our temperature results are equally ap- 
plicable to polariton arrays, making circuit- QED many- 
body systems equally difficult to realize.) The true na- 
ture of temperature in the phoniton array system and 
the potential for active cooling are subjects worthy of 
further consideration. Our proposed many-body systems 
with phonons can be developed further for the pursuit of 
quantum simulators [35 , 36 or mediators between differ- 
ent quantum components and potentially for new quan- 
tum devices. 
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